Interplay between Kondo suppression and Lifshitz transitions 
in YbRh2Si2 at high magnetic fields 
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We investigate the magnetic field dependent thermopower, thermal conductivity, resistivity and 
Hall effect in the heavy fermion metal YbRh2Si2. In contrast to reports on thermodynamic mea- 
surements, we find in total three transitions at high fields, rather than a single one at 10 T. Using the 
Mott formula together with renormalized band calculations, we identify Lifshitz transitions as their 
origin. The predictions of the calculations show that all experimental results rely on an interplay of 
a smooth suppression of the Kondo effect and the spin splitting of the flat hybridized bands. 



Since their discovery almost four decades ago, Kondo 
lattice systems pose a challenge to condensed-matter 
physicists. Above the characteristic Kondo temperature 
Tk they contain a periodic lattice of local 4/-derived 
paramagnetic moments. Below Tk, these moments be- 
come reduced and eventually fully screened by the con- 
duction electrons: The entanglement of the localized 4/ 
states with the delocalized conduction-band states leads 
to the formation of local Kondo singlets, which develop 
weak dispersion as a consequence of their periodic ar- 
rangement (Bloch's theorem) [1]. The delocalized Kondo 
singlets act as (composite) charge carriers, wihch exhibit 
a large effective mass (heavy fermions) due to the ex- 
tremely strong on-site Coulomb correlations. This is in- 
ferred from a huge Sommerfeld coefficient of the elec- 
tronic specific heat j2j but is invisible, e.g., in photo- 
electron spectroscopy [3], probing one-electron proper- 
ties. The localized 4/ states at T > Tr- are not part 
of the "small" Fermi surface (FS). However, being con- 
stituent to the composite charge carriers, the delocalized 
4/ states do contribute to a "large" FS below Tk- 

Recently, the effect of a magnetic field on these com- 
posite fermions has become an important issue. On the 
one hand most of the information on these quasiparti- 
cles comes from de Haas van Alphen (dHvA) experiments 
performed at high fields. However, band structure calcu- 
lations, necessary to analyze these experimental results, 
exist almost exclusively at zero field [H [5] . On the other 
hand, many heavy fermion systems undergo a metamag- 
netic transition at fields of the order of 10 T [BHS], where 
in some cases direct evidence (from dHvA) for a change 
of the FS and thus, of the heavy quasiparticles has been 
found [nilin], the origin of which remains controversial. 

YbRh2Si2 is such a case: Specific heat, susceptibility 
and magnetostriction measurements revealed anomalies 
at a critical field Bq = 10 T, indicating a rapid decrease 
of the effective mass at IS HU [E]. This was ini- 
tially interpreted as a breakdown of the Kondo screening 



resulting in a small FS on the high field side. Later, 
dHvA experiments were interpreted in terms of a Lif- 
shitz transition at i?o, where a spin-split band disappears 
[SI [13] . While these two interpretations seem quite differ- 
ent at first sight, they essentially rely on a similar model: 
Only in the presence of the very fiat bands of the com- 
posite fermions, the Zeeman splitting can induce such 
large effects on the FS at moderate fields as observed 
in YbRh2Si2. Simultaneously a splitting also drives the 
"large" FS continously toward the "small" one and leads 
to a decrease of the quasiparticle mass. Therefore, the 
field scale for both processes is related to Tk- 

In order to get more insight into the field evolution of 
the quasiparticles and into the connection between de- 
struction of the composite fermions and the occurrence 
of Lifshitz transitions, we performed detailed, field de- 
pendent transport studies in YbRh2Si2 in fields to at 
least 12 T. We focus on the effects far above the quan- 
tum critical point (QCP) &t B = 0.06 T which shall not 
be addressed here. In contrast to reports on thermody- 
namic measurements [9l [TTJ [12] , we show that the transi- 
tion at Bq consists of two close-lying features. Moreover, 
we observe an additional transition at 3.4 T. From the 
remarkable agreement between thermopower, electrical 
conductivity and magnetostriction above 2 T within the 
framework of Mott's formula, we demonstrate that all 
three transitions are caused by band structure effects. 

We compare our experimental results with predictions 
of the renormalized band (RB) structure calculation. 
They reveal three successive Lifshitz transitions at the 
three experimentally observed transition fields. This 
agreement demonstrates that the observed field evolu- 
tion results from an interplay of a smooth suppression 
of the Kondo effect and the spin splitting of a sharply 
structured density of states (DOS) generated by a strong 
anisotropic hybridization. 

We investigated three high quality samples from the 
same batch, with a residual resistivity of approximately 
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FIG. 1. Magnetoresistance, p, thermopower, S/T, and Lorenz 
ratio, L/Lo plotted as a function of magnetic field B. Lines 
and dots denote field sweeps, open squares are extracted from 
temperature sweeps, (a) p shows three kinks at B\ = 3.4 T, 
B2 ~ 9.3 T and B3 = 11.0 T. Their position is sample- and 
temperature-independent. Inset: enlarged view for 0.03 K. 

(b) S/T of sample ^1 at the same temperatures as in (a). 
The complex behaviour with several zero crossings reflects 
the multi-band character of YbRh2Si2. All three transitions 
are characterised by pronounced steps (insets: enlarged view). 

(c) The Lorenz ratio at 0.49 K is below unity with a minimum 
at about 0.5 T. 



0.5 fiflcm. On sample #1 we performed simultaneous 
d.c. resistivity, thermal conductivity and thermopower 
measurements down to O.IK and for 0.2T < B < 12T 
using the same soldered contacts, as described in Ref. HD 
Sample #2 was used to extend the resistivity data down 
to 0.03 K on another set-up applying an a.c. technique. 
We also performed Hall-effect measurements on sample 
#3 down to 0.05 K and up to 15 T. The magnetic field B 
was always applied perpendicular to the c-axis. The cur- 
rents for resistivity and thermal transport were parallel 
to B, for the Hall-cfFect perpendicular to c and B. 

Figures [ija) and [l|b) show the low-temperature mag- 
netoresistivity p{B) and the isothermal, field-dependent 



thermopower S{B), respectively. The Lorenz ratio 
L{B) /Lq = k{B)p{B)/TLo obtained from magnetore- 
sistance and thermal conductivity k using Sommerfeld's 
constant Lq = 7r^fc^/3e^, is shown in Fig. [ijc) only 
for 490 mK, because of an enhanced noise level at lower 
temperatures. Considering the B-T phase diagram of 
YbRh2Si2 m, it is natural to relate the low- field {B < 
2 T) behaviour of all these quantities to the signatures of 
the QCP and the surrounding non-Fermi liquid regime. 
Since our field range for the measurements on sample #1 
is limited to above 0.2 T, we observe only the high-field 
branch of these signatures, which extends to roughly 2 T. 
Nevertheless, our measurements agree with previously re- 
ported results: the step in magnetoresistance [15], the 
pronounced minimum in S{B)/T [16] and the minimum 
of the isothermal Lorenz ratio [T7] are clearly visible. 

Next, we focus on the high-field properties beyond 
2T, where YbRh2Si2 is a Fermi liquid below 200 mK 
[S]. The key features are three transitions visible as tiny 
kinks in the magnetoresistance at Bi — (3.4 ± 0.1) T, 
B2 ^ (9.3 ± 0.1) T and B3 ^ (11.0 ± 0.2) T in Fig. [ija). 
The transitions are more obvious in the thermopower 
(Fig. [l]^b)): S{B)/T shows three pronounced steps, 
which become sharper at lower T. 

These observations are interesting from two perspec- 
tives. Firstly, previous measurements have not observed 
features at 3.4 T P [Hi [12]. And secondly, the transi- 
tion reported in magnetisation, specific heat and mag- 
netostriction [SJ [TTl [T2] at roughly 10 T is actually com- 
posed of two well-separated features at 9.3 T and 11.0 T 
down to 30 mK, which do not merge in the extrapolation 

We further analyse our data using the Mott formula 
[l8] for the diffusion thermopower to clarify if the ob- 
served features have a thermodynamic origin. Phonon 
drag contributions to the thermopower are negligible in 
this temperature range [B] . The Mott formula generally 
holds at low temperatures in the absence of inelastic scat- 
tering [20j . To expand it, we first exchange the energy 
derivative of Incr with the field derivative using dB/de. 
In the second step, we use a = ne^r jm* together with 
m* cx (in 3D) and obtain 
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where a = l/p is the electrical conductivity, n the to- 
tal electron concentration, N the DOS, m* the effective 
mass, and r the scattering time. This splits the ther- 
mopower into a scattering part (r) and a part represent- 
ing the band structure (TV). 

To test the first expansion in Equation [l} we compare 
S{B)/T with dlna/dB in Fig. [ga). Above 3T, both 
curves can be scaled on top of each other using a con- 
stant dB/de. This suggests a linear relationship between 
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e and which we write in terms of an effective Zee- 
man energy e = geg^sB /2, with ^eff of 16 ± 1 (hb is 
the Bohr magneton), and discuss later on. Moreover, it 
implies that inelastic scattering (which would invalidate 
the Mott formula) is insignificant in this regime. In the 
low-field region below 2T both curves show a disparate 
behaviour. Consequently either e oc i? is violated or in- 
elastic scattering is significant. Both are likely to occur 
close to a QCP: The first can arise if there are dramatic 
changes in the magnetisation, the band structure or the 
FS. Inelastic scattering can increase close to a phase tran- 
sition. 

We now turn to the second expansion of Equation [l] 
and use the linear magnetostriction coefficient A as a 
thermodynamic probe to access the contribution from 
the DOS, since the data available have a higher res- 
olution than e.g. specific heat or magnetisation. The 
inset in Fig. ^h) demonstrates that the linear magne- 
tostriction {Al \\ B ± c) follows Al/l (X AP-'^ over the 
entire field range. In a Fermi liquid x — dM/dB fol- 
lows the behaviour of the DOS, and therefore also does 
A = d{Al/l)/dB. Applying a constant scaling factor, 
dlnX/dB matches S/T in Fig. ^h), even at low fields. 
Only the double hump around 10 T is more pronounced. 
Considering the nice agreement at low fields, it is likely 
that the discrepancy of S and a in the same low-field 
regime is due to enhanced inelastic scattering rather than 
a failure of e ~ B. The surprisingly good qualitative 
agreement, especially at Bi, B2 and B3, is clear evidence 
that the origin of all three transitions lies within the cor- 
related band structure of YbRh2Si2. 

Our experimental data together with previous results 
already indicate the complex nature of these band struc- 
ture effects. For example, the transition around Bq is 
composed of two fields and the drop in the specific heat 
coefficient ^{B) between them is only moderate with a 
decrease from about 250mJ/molK^ to lOOmJ/molK^. It 
is therefore unlikely that either a single Lifshitz transition 
[S] or a sudden localization of the /-electrons [21 HH [12] 
is — on its own — a sufficient model to describe the be- 
haviour around 10 T. Hence, a theory of the field evo- 
lution of YbRh2Si2 should include not only the Kondo 
effect to describe de-renormalization processes, but also 
the specific correlated band structure of YbRh2Si2 to re- 
veal topological transitions. 

We therefore conducted the field-dependent RB cal- 
culations summarized in Fig. [3| which are an extension 
to the results from Ref. [H] with a focus on the detailed 
development of the FS and the DOS in field. We used 
the RB method described in Ref. (211 and [H| with a very 
tight fc-mesh in zero field of 8125 points in the irreducible 
wedge to resolve changes in the isoenergy surfaces. In fi- 
nite fields we used 405 /c-points. 

Figure [Sj^a) shows the calculated zero-field quasipar- 
ticle DOS N{e) and Fig. [sjjb) displays the correspond- 
ing isoenergy surfaces of the most important FS sheets 
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FIG. 2. Analysis of the thermopower with the Mott for- 
mula. The measured SjT (dots) is compared with that cal- 
culated (lines) from (a) the electrical conductivity and from 
(b) the magnetostriction coefficient using Equation [T] with 
d\nN/dB cx d\n\/dB and dlnr/dB = 0. Within the Fermi 
liquid regime (B > 2 — 3 T) all curves show the same overall 
behaviour. Inset: Magnetisation dependence of the magne- 
tostriction. The data (solid blue) perfectly follows a power 
law Al/l oc M^'^ (dashed black line). 



with /-character — the so called "doughnut" (top) and the 
"jungle gym" (bottom), respectively. Scanning through 
N[e) within the displayed energy range, the calculations 
predict the four color-coded regimes characterized by dif- 
ferent topologies of the isoenergy surfaces, separated by 
Lifshitz transitions. 

An applied magnetic field spin splits the DOS into 
a minority and a majority branch. These do not shift 
rigidly in field but the amplitude of the coherence peak 
decreases (Fig.jsjjc)), mainly because of the weakening of 
the Kondo effect. Nevertheless the characteristics of the 
band structure are not changed. Therefore we expect the 
energy evolution of the isoenergy surfaces in Fig. [sjjb) to 
be the same as the field evolution of the FS. I.e., the FS 
of the majority band stays in the red (D) regime, while 
the FS topology of the minority band changes from the 
red (D) type through yellow (C) and green (B) to the one 
marked in blue (A). 

We take advantage of the similar shape of the zero- 
and finite-field DOS and assign the energy of a topology 
change (in zero field) to a magnetic field where the corre- 
sponding feature in the DOS crosses ep: The transition 
from the blue (A) to the green (B) regime corresponds to 
a kink in the DOS which reaches ep for B^ = (11 ± 1) T. 
Similarly, we obtain B2 — (9± 1) T for the second transi- 
tion (green (B) - yellow (C)). The difference B3 — B2 fits 
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FIG. 3. Isoenergy surfaces for 5 = and the quasiparti- 
cle DOS development in finite field calculated using the RB 
method. The zero-field DOS in (a) is divided into four re- 
gions (blue (A), green (B), yellow (C), red (D)) distinguished 
by different topologies of the isoenergy surfaces shown in (b). 
In the yellow region (C), we show in plan the "jungle gym" 
FS sheet exactly at the topological transition between B and 
C. (c) illustrates the magnetic field evolution of the DOS for 
selected fields, with the zero-field DOS in grey for compari- 
son. Inset in (a): Energy-field map of the DOS interpolated 
in 1 T steps. One can assign the four energy regions and their 
isoenergy surfaces in a) and b) to the four field ranges and 
their FS separated by Bi, B2 and Bs. 



to the linear shift of O.lmeV/T (inset Fig. [3 a)), which 
in turn matches the ESR 51- factor of 3.5 [23 (applying 
e — gfiBB/2). We use this shift to estimate the field cor- 
responding to the weak third transition from yellow (C) 
to red (D) to = (2.5 ± 1) T. These fields, extracted 
from the calculations, are in excellent agreement with the 
transitions B2 and B3 found in the experiment. This 
proves that they correspond to three Lifshitz transitions 
of the types illustrated in Fig. [sjjb) . 

Additionally, the linear sweep of the van-Hove singu- 
larity predicted by our calculations confirms the linear 
energy-field dependence found in our data analysis. It 
is worth reiterating that the only adjustable parameter 
leading to the remarkable accuracy of the Mott formula 
in the two comparisons shown in Fig. [2] is gcs- Moreover, 
the thermopower as well as d\nX/dB and dlna/dB are 
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FIG. 4. Hall effect on sample #3. The field-dependent inverse 
Hall coefficient shows several extrema close to Bi,B2 and B2. 
obtained from resistivity data (see Fig.[T][a)). Inset: the Hall 
resistivity, pxy, is almost linear below 0.2 K. 



independent of sample geometry, thus systematic errors 
are almost negligible. The difference between goff = 16 
(corresponding to a de/dB ~ 0.5meV/T) and g — 3.5 
(O.lmeV/T) from our calculations, however, is unsur- 
prising, since a rigid band shift is obviously insufficient to 
account for the experimental data. Unexpectedly, these 
field-induced changes of the band structure also enter 
linearly into gcff, a fact which certainly deserves further 
studies. 

Our results confirm on the one hand that the Fermi liq- 
uid state observed in YbRh2Si2 in finite magnetic fields 
originates from the formation of local Kondo singlets. A 
magnetic field supresses the local Kondo effect. How- 
ever, Numerical Renormalization Group method for a 
single impurity (see e.g. [MUTT] ) predicts a continous de- 
crease of the effective mass. Therefore, our calculations 
demonstrate on the other hand that the single impurity 
model alone cannot account for the abrupt changes in 
the observed quantities reported here. The latter have 
to be attributed to coherence effects arising from the pe- 
riodic arrangement of the Kondo ions. The anisotropic 
hybridization of the 4/ states with the conduction bands, 
caused by the highly anisotropic crystalline electric field 
(CEF) ground state, leads to van-Hove-type singularities 
in the quasiparticle DOS. The structures in the quasi- 
particle DOS highlight changes of the iso-energy surfaces 
as shown in Fig. |3] In a magnetic field, the 4/ states 
are split which, in turn, leads to a Zeeman splitting of 
the quasiparticle bands. The relative shifts of the latter, 
however, are enhanced by a field-dependent Sommerfeld- 
Wilson ratio which reflects the local many-body effects. 

Experimentally, further insight into the evolution of 
the FS can be obtained by Hall effect measurements. The 
isothermal Hall resistivity Pxy{B) for selected T is pre- 
sented in the inset of Fig. |4] and indicates hole- dominated 
transport. Importantly, beyond 12 T all curves coincide 
pointing towards a field-driven suppression of the local 
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Kondo effect. This suppression is a continuous process 
and likely accounts for the maximum in the effective car- 
rier concentration, l/i?H ~ V^^zg^' around 6 T and 
lowest temperatures via the hybridization of conduction 
and 4/ electrons. For T > 0.2 K an anomalous contribu- 
tion [55] to i?H comes into play (c/. also [53]) such that 
1/Rh no longer tracks the carrier concentration. Quali- 
tatively, the transitions at Bi , B2 and B3 are observed as 
extrema in Ryi{B) although — quantitatively — at slightly 
different fields. The latter might be caused by the com- 
plex nature of the Hall effect [50] {e.g. multiple bands) or 
the different current alignment j -L B {j\\B for all other 
measurements) . 

In conclusion, we showed that the thermopower is an 
excellent tool to reveal field induced changes in the FS 
of a correlated metal, hardly detectable by any other 
probe. In YbRh2Si2, we find three successive transitions, 
which were identified as Lifshitz transitions by a com- 
parison with predictions from detailed RB calculations. 
This implies, that the unusual high-field properties of 
YbRh2Si2 arise from the interplay of (a) the symmetry 
of the CEF ground state (g-factor, anisotropic hybridiza- 
tion), (b) the suppression of the local Kondo effect (re- 
duced effective mass, field-dependent Sommerfeld- Wilson 
ratio) and (c) the coherence effects due to the periodic- 
ity of the lattice (van-Hove singularities, Lifshitz transi- 
tions). The excellent agreement between our experimen- 
tal results and our RB calculations demonstrates that 
RB calculations are a very suitable approach to describe 
quasiparticles in the Kondo lattice. 
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